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Abstract. 

We study the robustness of the spectrum emitted by an acoustic black hole by considering 
series of stationary flows that become either subsonic or supersonic, i.e. when the horizon 
disappears. We work with the superluminal Bogoliubov dispersion of Bose-Einstein condensates. 
■ We find that the spectrum remains remarkably Planckian until the horizon disappears. When 

O" 1 the flow is everywhere supersonic, new pair creation channels open. This will be the subject of 

| , a forthcoming work. 

> . 

, 1. Introduction 

' In the hydrodynamic approximation, i.e. for long wavelengths, the propagation of sound waves 

in a moving fuid is analogous to that of light in a curved spacetime [H [2] . More precisely, a one 
dimensional flow defines the acoustic metric 

(N 

ds 2 = -c 2 dt 2 + (dx-vdt) 2 , (1) 

where v is the flow velocity and c the speed of sound. Assuming that the fluid flows from right 
to left (v < 0), a sonic horizon is present where w = c + v crosses 0. This situation describes 
rN ■ a black (white) hole horizon when kk = d x w evaluated at w = is positive (negative). One 

thus expects that such system would emit a Hawking flux, i.e. a thermal flux of phonons at a 
temperature given by Th = kr/27t (in units where ks = h= 1). 

However, since Hawking radiation relies on short wavelength modes [3J 0], the dispersion 
of sound waves, which is neglected in the hydrodynamic approximation, must be taken into 
account. Following the original work [S], this was done using analytical El and 
numerical \10\ \TT\ [T2"| [T3] methods. Provided the dispersive scale £, the "healing length", is 
much smaller than the surface gravity scale c/n, it was found that the thermal properties of the 
spectrum are extremely robust. In fact, so robust that it is very difficult to characterize and 
analyze the small deviations with respect to the standard flux. 

Our aim is to complete the analysis of |12t [T3] where it was shown that the spectrum remains 
Planckian whenever Th < cj max /15, where w max is a critical frequency which scales as l/£ but, 
more importantly, also depends on w(x = — oo), the asymptotic value of c + v. Since this result 
was obtained using symmetrical profiles like w ~ D t&nh(Kx/D), one may wonder whether the 
robustness of the spectrum was not partially due to the symmetry of the flow w.r.t. the horizon, 
namely w(—x) = —w(x). To investigate this question, we shall consider flows which are highly 
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Figure 1 . Four velocity profiles w = c + v of Eq. ([2]) with the same height D but four values 
of D\. For D\ = (solid line), one recovers the symmetric case of |12j with the Killing horizon 
(solid vertical line) at x = 0. For —D<D\ < (dashed line), the horizon (dashed vertical line) 
is shifted to x > 0. In the last two cases, there is no horizon: for D\ < —D (dotted line), the 
flow is everywhere supersonic, and for D± > D (dotdashed line), it is subsonic. 

asymmetric. In particular we shall study the spectrum in the interesting limiting cases where 
the horizon disappears because the flux becomes either sub or supersonic. In parallel with this 
work, we investigated in p3] other classes of profiles generalizing those of |12|. [13] . 

2. Phonon spectra in asymmetric flows 

We consider elongated condensates, stationary flowing along the longitudinal direction x. We 
assume that the transverse dimensions are small enough that the relevant phonon excitations 
are longitudinal. The system is thus effectively described by the 1 + 1 dimension metric (UJ 
where v and c depend only on x. We work with flows w = v + c given by 



D fixes the gap [io(oo) — w(— oo)]/2co, while D±, which was taken to be in [12] [13], fixes 
the asymmetry between sub and supersonic regimes. When D\ > D, the flow is everywhere 
subsonic, and when D\ < —D, it is supersonic, see Fig. [TJ Instead, when \D\\ < D, there is a 
sonic (Killing) horizon where w = 0, localized at 



It separates the sub (x > xh) from the supersonic region (x < xh). In such a flow, when ignoring 
dispersion, the spectrum of upstream phonons spontaneously emitted from the horizon would 
be very simple, and would strictly correspond to the Hawking radiation [TJ [5]: It would follow 
a Planck law at Hawking temperature Th = kk/27t. The surface gravity kk is here 



Taking dispersion into account, the spectrum becomes much more complicated. In this paper, 
we consider phonons in Bose-Einstein condensates (BEC) with the Bogoliubov dispersion 
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Figure 2. Left panel: T u for 6 values of D\/D, with D = 0.5 and k£/\/2c = 0.1 fixed. The 
horizontal lines represent Tn(D\) of Eq. respectively for \D±\ = 0,0.2,0.4. For D\ < —D the 
curve is truncated at w m i n as explained in the text. Right panel: Ah as a function of w max /K for 
two values of D\/D and various D\. When w max /ft is large enough, deviations from Planckianity 
are small and, more importantly, governed by the ratio D\/D only. 



where uj is the lab conserved frequency. In the sequel, following the numerical procedure of [13] . 
we study the deviations of n u , the phonon spectrum, w.r.t. the standard Planck one, by varying 
the offset parameter D\ while keeping fix all the other parameters. 

To characterize these deviations, we use the temperature function T u defined by 

n u> = , L v t, (6) 

exp(w/r w ) - 1 

which is constant when the spectrum is Planckian. This quantity thus allows a direct comparison 
with the Hawking temperature Th = k^/2tt of Eq. In Fig. [21 T w and Th (horizontal lines) 
are plotted for various values of D\ . We distinguish four regimes: 

• For \D\\ sufficiently smaller than D, the spectrum is Planckian until uj approaches w max , 
where it vanishes. For D\ > 0, w max decreases since \w(x = — oo)| does so, and conversely 
for D\ < 0. However for both signs, at low frequency T u closely follows kk/2vt of 

• When \D\\ —> D. kk of drops down to and the sonic horizon disappears. In this 
critical regime, deviations from Planckianity appears even at low frequency. 

• For D\ > D, the flow is everywhere subsonic. There is no particle production because there 
are no negative norm modes with uj > 0. 

• For D\ < —D, the flow is supersonic. A new critical frequency w m ; n < w max appears. For 
< uj < w m ; n , there are now 4 asymptotic in and out modes, and the scattering matrix 
is 4 x 4. Since the code of [T3] is designed to the 3x3 case, this regime requires a new 
code [TS]. However, for oj m \ a < uj < a/ max there are only 3 modes, and our code can handle 
this frequency band. In this range, as can be seen in the left panel, is close to the case 
where D\ > —D, even though the spectrum is not at all Planckian. 

Furthermore, we checked that the shape of the spectrum is hardly changed when varying the 
healing length £. Namely, spectra calculated with different values of £ almost coincide when 
T u is plotted versus the rescaled quantity w/w max (£). To identify the parameters governing the 
deviations from Planckianity we define 

A _ /(^ = t h) ~ /t h (^ = ?h) , , 

AH = /*> = Th) ' (7) 



where f(ui) = un u is the actual energy flux and /t h ( w ) i s the thermal flux at temperature Th- 
In Fig. [21 right panel, Ah is plotted as a function of w m ax(£)/ K - Various values of D\ are used, 
while the ratio D\/D is kept constant respectively for the three lower curves {D\/D = —1/5), 
and for the three upper curves {D\/D = 4/5). As expected, Ah goes to for large values of 
Wubx/k. What we learn here is that the leading deviations from Planckianity do not depend 
separately on D\ or D but only on their ratio. 

3. Comments and conclusions 

In this work, we have studied the robustness of the properties of black hole radiation in BEC 
by considering flows characterized by a high asymmetry between the sub and the supersonic 
region. We found that the spectrum remains Planckian until one approaches the critical cases 
(|-Di| — > D) where the surface gravity vanishes. 

In fact, for < D\ — > D, since the fluid flows in the supersonic region with a velocity just 
above that of sound, the cutoff frequency uj max goes to zero. Hence the range of frequencies 
where the temperature T u is nearly constant shrinks and disappears in the limit. As expected, 
no emission of particles is found when the flow remains everywhere subsonic. 

The situation is much more interesting in the opposite case, for > D± — > —D, when pushing 
the velocity of the subsonic region towards the speed of sound. In this regime, even though kk 
of (jl]) goes again to 0, u; max now increases. As a result, the low frequency plateau of T w at 
low temperature again shrinks, but the spectrum becomes more and more blue. This can be 
understood from the fact that more dispersive modes, i.e. modes with higher cj, have still their 
turning point near x = where the gradient of w is K S> kk- 

When D\ < —D, there is no sonic horizon since the flow is everywhere supersonic. Yet 
the spectrum is continuously deformed even though the Planck character is completely lost in 
the frequency range w m i n < oj < w max , where oj m \ D is a new critical frequency. For u < w m m, 
new scattering channels open because there is a fourth asymptotic mode. Investigations of this 
case [J5] are in progress. 

To conclude, we stress that the thermal spectrum at the naive surface gravity ([3]) provides 
a reliable approximation of the actual phonon spectrum in all cases but the critical ones where 
the sonic horizon disappears. In these extreme cases, the spectrum is no longer thermal and its 
properties are governed by the dispersive properties of field. 
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